Stokes— Einstein relation violation 
and the upper critical dimension of the glass transition 
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Using numerical simulations we investigate the dimension evolution of the Stokes-Einstein relation 
violation in a simple glass forming system. The results strongly support the existence of an upper 
critical dimension d u = 8 above which the Stokes-Einstein relation is not violated, as predicted by 
the mean-field theory of glasses. In addition, the associated critical exponents evolve linearly in 8 — d 
below d = 8, which suggests that a suitable renormalization group expansion around the mean-field 
theory could provide a quantitative description of the glass problem. An analysis of single-particle 
motions in terms of Levy flights offers a possible scenario for relating these corrections to the mean- 
field description with the dynamical heterogeneity of glass formers. 
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The glass problem is notoriously difficult. Construct- 
ing a completely satisfying theory for how rapidly a fluid 
turns sluggish with only unremarkable structural changes 
remains hotly debated [1] . A key hurdle is that the theo- 
ries of the glass transition typically provide insufficiently 
precise predictions for decisive experimental or numeri- 
cal tests to be performed. Identifying the "smoking gun" 
that could tell theories apart has thus long been a prior- 
ity. 

A promising way of addressing this issue is to investi- 
gate the glass transition as a function of spatial dimen- 
sion d, which was recently achieved in numerical simula- 
tions [2] . In low d various phenomena kinctically compete 
in glass- forming fluids: (i) crystal nucleation, (ii) barrier 
hopping due to thermal activation, and (Hi) trapping in 
phase space due to the proximity of crgodicity breaking. 
Different theories give more or less weight to these phys- 
ical processes, which are so well enmeshed that they are 
typically hard to tell apart. An advantage of increasing d 
is that both nucleation and barrier hopping get strongly 
suppressed. Attention can then be focused on the onset 
of ergodicity breaking [2] . 

A central quantity on which most theories give spe- 
cific predictions is the violation of the Stokes-Einstein 
relation (SER). According to the SER the structural re- 
laxation time T a , which is proportional to the viscosity 
in a dense fluid, and diffusivity D should give a constant 
Dr a . Yet in d=2 and 3, a large violation of this relation 
is observed upon approaching the glass transition [3-6]. 
In most cases it is found that D oc t~ 1+u [7-9]. The SER 
violation is traditionally attributed to a large spatial het- 
erogeneity of the "local relaxation times" [9]. Diffusion 
is dominated by the faster regions, while viscosity and r Q 
are dominated by the slower regions. Hence, 1/D <C r a 
and one expects lu > [9] , as is indeed observed numeri- 
cally [6-8]. 

Explanations for the evolution of the SER violation 
with d fall into two categories. First, theories based on 
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FIG. 1: Mam panel: The self-diffusion D and the structural 
relaxation time r a are related by a power-law with exponent 
— 1 + ui (dashed lines), where w changes with d. Note that the 
power-law regime is truncated at low d by activated events on 
a scale similar to that observed for other dynamical quanti- 
ties [2]. Over the limited remaining dynamical rang for d=3, 
tj=0.19(3). Inset: The exponent u (symbol) linearly vanishes 
as u = A(d u - d) with d u = 8 and A = 0.038(3) (solid line). 



adapting the standard description of critical phenomena 
to the glass problem predict an upper critical dimension 
d u - For d > d u a mean-field theory is expected to give 
the correct critical description, while for d < d u criti- 
cal fluctuations qualitatively change the behavior of the 
system by renormalizing the scaling exponents or elim- 
inating the transition altogether. In this scenario, SER 
violations would disappear for d > d u , and it is expected 
that d u = 8 based both on static [10, 11] and dynam- 
ical [9] descriptions. In particular, a scaling argument 
proposed in Ref. [9] gives uj ~ (8 — d)/(4'y), where 7 is 
the exponent that describes the divergence of the relax- 
ation time with density (and, more generally, tempera- 
ture) near the onset of ergodicity breaking [2]. Second, 
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FIG. 2: Main panel: As the system becomes more sluggish. 
P{r,T a ) (here for d=4) develops a long tail that is partly 
hidden by structural oscillations on a scale r ~ 1. Inset: 
As dimension increases, the tail of P(r,T a ) for isodiffusive 
systems (2dD = 0.0008(3)) gradually disappears. 



other theories, e.g., those based on a certain class of ki- 
netic models [12], predict that the critical behavior of 
the system is qualitatively similar in all spatial dimen- 
sions, hence d u — oo [7]. In the East model studied in 
Refs. [7, 12, 13], for instance, it is found that ui « 0.2 for 
all d. The existence of d u can thus be directly established 
by studying the SER violation as a function of d. An at- 
tempt in this direction was done in Refs. [8, 14]. In this 
letter, we perform numerical simulations to measure ui in 
higher dimensions, and find that d u = 8. We also ratio- 
nalize that the SER violation can be interpreted in terms 
of a Levy flight that describes the underlying behavior of 
the particle displacement distribution on the scale of r Q . 

Simulations - Event-driven molecular dynamics sim- 
ulations of hard spheres in dimension d=3-9 are per- 
formed, following the method of Refs. [2, 15, 16] for 
at least N = 8000 identical particles of unit diameter 
<7, except for d=3 where an equimolar binary mixture 
with a diameter ratio 7:5 is used to prevent crystalliza- 
tion [17, 18]. Simulations are run at constant unit inverse 
temperature j3 with time t expressed in units of \J (3ma 2 
for particles of unit mass to. The diffusivity is extracted 
from the long-time behavior of the mean-square displace- 
ment D = limt^ oc (r 2 ) t /2(ii, where r = r(t) — r(0), and 
r a is the characteristic decay time of the self-intermediate 
scattering function F s (k*,t) = (e lk *' r ) at the wave vec- 
tor k* that maximizes the first peak of the structure fac- 
tor. Figure 1 confirms that upon approaching the glass 
transition the regime where D oc t~ 1+u extends with 
dimension, in agreement with activated relaxation pro- 
cesses getting increasingly suppressed. The results are 
also consistent with u vanishing linearly at d u = 8, which 
indicates that a finite upper critical dimension indeed ex- 
ists. However, we find u> ~ 0.038(8— d), while the analysis 
of Rcf. [9] gives uj ~ 0.12(8 — d), using earlier evaluations 



of 7 « 2.1 [2]. This discrepancy calls for a better inves- 
tigation of the physical origin of the SER violation for 
d < d u - 

Both D and r Q arc themselves functionals of the self 
van Hove function G s (r,t) [19]. Because D weighs large 
displacements more than r Q , the SER violation thus sug- 
gests that G s {r,T a ) has an anomalous behavior at large 
r. The spherically integrated displacement distribution 
P(r,T a ) = Qdr d ~ 1 G s (r,T a ) indeed develops a fat tail at 
large r when t q increases (Fig. 2), as was previously re- 
ported in various d=2 and 3 systems [18, 20, 21]. This 
behavior is remarkably distinct from the Gaussian dis- 
tribution expected in the long-time diffusive regime [19], 
and is also different from the exponential regime that is 
observed in the caging, or f3, regime [2]. In the struc- 
tural relaxation regime, the fat tail of G s (r, r Q ) results 
in a larger than expected mean square displacement and 
in the SER violation. By contrast, with increasing d the 
bulging tail resorbs and the SER is gradually recovered 
(Fig. 2). To further substantiate this claim, we introduce 
a time scale td, such that (r 2 ) TD = 2d/(k*) 2 . Numeri- 
cally, we find that (i) td corresponds to the beginning of 
the diffusive regime and thus (r 2 ) t ~ 2dDt for t > Tr>, 
and (ii) D oc r^ 1 (Fig. 3). The SER violation and (ii) 
imply that rrj ~ Tq~" *C T a . Hence from (i) wc deduce 
(r ) Ta ~ 2dDr a oc T a /tjj oc t£, i.e. the second moment 
of G s (r, r a ) diverges with r a . This argument establishes 
the relation between the SER violation and the fat tail 
of the self van Hove function, as illustrated in Fig. 3. 

Levy Flight Model - A minimal model for describing 
the physical consequences of this distribution anomaly 
can be formulated in terms of a truncated Levy flight [22- 
24], which amount to the following two assumptions, (a) 
For large r Q , the distribution P(r, r Q ) converges to a finite 
limit around its maximum ?*t yp (T Q ) ~ 1/fc* that is almost 
independent of r a (Fig. 2), while its behavior for r 
r typ( T «) i s described by 

P(r,r a ) = r -^^T[r/ar a )}, (1) 

where we defined a dynamical length £ ~ r« and the 
critical exponents /j, and C > 0. In this model "P(O) is 
a constant and V{x) rapidly vanishes at large x. We 
assume that fj, < 2 for d < 8, in order for the second mo- 
ment of P(r, T a ) to diverge when r a — > oo. Note that the 
bounded limit distribution at small r guarantees its nor- 
malizability, which is consistent with the majority of dis- 
placements taking place on a scale rt yp (T Q ) at r Q . (b) On 
a time scale f>T„, particle displacements are described 
by a random walk with uncorrelated elementary steps 
distributed according to P(r, r a ). Stated differently, the 
self van Hove function at r = tjr a > 1 can be obtained 
by convolutions of the function evaluated at r = 1, and 
thus 

F s (M)~e-^«( fc2 ) . (2) 
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FIG. 3: Upper panel: The mean square displacement (r 2 )t = 
R 2 (2,t) as a function of time for d — 5 and several densi- 
ties upon approaching the glass transition. Crosses indicate 
the time td such that (r 2 ) TD = 2d/(k*) 2 (long-dashed line), 
which marks the beginning of the diffusive regime. Asterisks 
indicate the relaxation time T a , which shows the growth of 



(short-dashed line). Inset: The separate evolu- 



tion of r a (asterisks) and to (crosses) with D, where D ~ r £ 



(long-dashed line) and D 



(short-dashed line). Lower 



panel: R (z,t) for z=0.1, 2, and 4 in d=5 for ip = 0.268 
(squares). The Levy flight regime extends from t > r„ to the 
rightmost part of the figure, where the functions cross over 
to the Gaussian diffusive regime. The predictions of Eq. (7) 
are consistent with the simulation results. The dot-dashed 
line has slope 2/fj, with fi — 1.423, following Eq. (6), while 
the full line has slope 1. The non-Gaussian parameter oi2(t) 
(circles) shows that non-Gaussianity is maximal at td and 
persists across the Levy flight regime, while aa(^) ~ 1/t at 
larger times in the Gaussian diffusive regime. 



If these two assumptions are true, at large times particles 
perform a random walk whose distribution of elementary 
steps acquires a power-law tail at the onset of ergodic- 
ity breaking, i.e. when the structural relaxation time 
diverges. This Levy flight description of the displace- 
ment tail at r a offers a natural microscopic explanation 
for SER violation. Numerically, however, the large r tail 
of P(r,r a ) is modulated by structural oscillations that 
partially obscure the presumed power-law behavior (see 
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FIG. 4: Main panel: The value of a2(r a ) versus r a for several 
dimensions indicates that 02 (r a ) ~ t2 over the non-activated 
regime. Inset: The exponent 77 (symbol) scales linearly 77 = 
A(d u -d)+jjoo with d u = 8, A = 0.017(2), and 7?oo = 0.23(1) 
(solid line), which suggests that 77 ~ 0.315 in d=3. 



Fig. 2). More robust consequences of the model must 
thus be examined to assess its relevance. 

Exponents and Scaling - Considering the fractional 
moments of the displacements, defined as R 2 (z,t) = 
[(r z )t]~ , one can show from Eq. (1) that 

RH? rlJ n yp (T a )^/^(r a ) 2{1 ^ /z) for z > /i , 
1 ' I r typ (r Q ) 2 C(z) forz< M) 

(3) 

where C{z) remains finite for r Q — >• 00 (but diverges for 
z — > fi). Using this form, and recalling that fj. < 2, we 
obtain that 

{r i) Ta = i? 2 (2,r Q ) ~ e(r Q ) 2(1 ^ /2) ~ r^ 1 -^ . (4) 

Hence, uj = 2(1— /i/2)/£andthe non-Gaussian parameter 
scales as 



d (r 4 ) Ta 
d+2<r 2 ) 2 



[R\^r a )Y /c 

[i? 2 (2,T Q )] 2 a 



(5) 

Numerical results give 77 - 0.23(1) + 0.017(2)(8 - d) 
(Fig. 4), so from the measured d < 8 exponents we de- 
duce that 



0.46 + 0.034(8 - d) 
+ ~ 0.23 + 0.055(8 - d) ' 
1 77 + uj 1 , 



(6) 



■[0.23 + 0.055(8 -d)} 



which gives /i sal. 25 and 1/C «0.25 in d=3. 

From Eq. (2) we also obtain that for finite r Q there 
exists a "Levy flight" regime 1 < t -C £ m i n which 

R*(z t)~ / Wr) 2 ^ 1 -^ forz>M, (?] 
{ , ) I n yP (r) 2 C(z) forz< /J , t?) 
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with a typical displacement rt yp (r) ~ t 1 ^. When r 



and thus t ~ r a ^ 



-1+V 



we obtain R (z,t) 



r typ (^) 2 ~ £ 2 for all z. From that point on, i.e. for 
t > t^ +v , the two length scales r ty p(r) and £ merge and 
one recovers the trivial Gaussian diffusive regime. The 
curves for all z then merge and scale as r. Note that 
Eq. (f) is valid for r Q < t < t^ +v , with a i-dependent 
r typ . Although our results are well captured by keeping 
£ constant, we cannot exclude the possibility of a small 
time dependence of £ in this regime. 

Figure 3 shows that R 2 (z,t) follows the prediction of 
Eq. (7) over one decade of t in d = 5, which strongly 
suggests that the postulated model is accurate. Be- 
cause we can only access a limited dynamical range and 
the exponent ?/ is quite small, the Levy flight regime 
r Q < t < T^ +n is, however, rather limited and therefore a 
complete verification of Eq. (7) is not yet possible. The 
key difficulty is that the range of accessible r Q is always 
limited. Activated events interfere with the slowdown in 
<i=3-5, and numerical simulations of denser systems in 
d > 6 are currently beyond computational reach. 

Conclusion - In summary, the critical behavior of 
the SER violation disappears for d > 8, which indicates 
that the glass transition has an upper critical dimension 
d u = 8. For d > d u the mean-field theory of the glass 
transition is thus a good description, and extending this 
theory using renormalization group methods should top 
the research agenda, as suggested in Ref. [9]. In par- 
ticular, a better understanding of the physical behavior 
below d u would be helpful to clarify why the prediction of 
Ref. [9] for the exponent uj is not quantitatively accurate. 

Of particular interest in this respect is understanding 
how dynamical heterogeneity, which is well documented 
in the literature [8, 18, 25] and is confirmed by a prelim- 
inary analysis of our data in all dimensions, results in a 
fat displacement tail at r Q . Our analysis of particle dis- 
placements only covers the regime t > r a , so we still lack 
a complete understanding of the microscopic mechanism 
on scales t < r Q . The question is especially interesting 
around tu, where diffusion sets in and non-Gaussianity is 
maximal. A possible explanation for the Levy flight-type 
distributions is for displacements on shorter timcscalcs 
to be random in time, but correlated in space, result- 
ing in large spatial jumps [24]. Because this process is 
necessarily collective, the description also offers a direct 
link between the SER violation and a growing dynamical 
heterogeneity [24]. 

A related point that deserves further examination is 
the connection between the dynamical length £ from the 
Levy distribution expression in Eq. (1) and the spatial 
extension of dynamical heterogeneity £4, as quantified 
by four-point correlations [8, 18, 25]. An argument pro- 
posed in Ref. [24] relates these two lengths in a differ- 
ent dynamical regime, but it might be possible to ex- 
tend the argument to the regime considered here. Mean- 
field treatments indeed predict that £4 ~ tU^ 1 ^ for d > 



8 [10, II, 26], which, using 7 = 2.1 [2], gives £4 ~ t°- 12 . 
This scaling is remarkably close to the one we obtain here 
for the dynamical length, i.e. £ ~ Ta = t^ 115 in d = 8. 
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